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Abstract. The Gelfond-Khovanskii residue formula computes the sum of the 
values of any Laurent polynomial over solutions of a system of Laurent poly- 
nomial equations whose Newton polytopes have sufficiently general relative 
position. We discuss two important consequences of this result: an explicit 
elimination algorithm for such systems and a new formula for the mixed vol- 
ume. The integer coefficients that appear in the Gelfond-Khovanskii residue 
formula are geometric invariants that depend only on combinatorics of the 
polytopes. We explain how to compute them explicitly. 



1. Introduction 

A Laurent polynomial is a finite linear combination of monomials with complex 
coefficients: 

f(t)=Y, Xktk > t k = l$ 1 ...tfr, k = (k 1 ,...,k n ), X k eC. 

The convex hull of those lattice points k G Z" for which ^ is a convex polytope 
in K", which is called the Newton polytope of /. The value of /(<) is defined for 
all t in the set (C \ 0)" = {(eC" \ h . . . t n ^ 0}. 

In this paper we talk about systems of n Laurent polynomial equations in 
(C \ 0)™ whose Newton polytopes have sufficiently general relative positions (the 
precise definition will be given in Section 2.1). O. Gelfond and A. Khovanskii 
[GK02] discovered that one can say a lot about solutions of such systems explicitly 
in terms of the Newton polytopes and the coefficients of the system. The following 
theorem ([GK02], Corollary 1.6) illustrates that. 

Theorem 1. Consider a system of n Laurent polynomial equations 

(1) = ••• = /« (*)=0, fe(C\0) n , 

whose Newton polytopes Ai, . . . , A n have sufficiently general relative positions. The 
sum of the values of a Laurent polynomial q over solutions 1 of (1) counted with 
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multiplicities is equal to 

(- 1 ) 2^ c A res A ^— A •••A— J, 

Aevort(A) •' 1 ■'" 

w/iere ffte sum is taken over the vertices A of the Minkowski sum A = Ai + - • - + A n , 
vesA(q^- A • • • A is the residue at the vertex A, and ca is the combinatorial 
coefficient at A. 

The residue at a vertex A is an explicit Laurent polynomial in the coefficients 
of q and /i, . . . , /„. In the case when n = 1 the polytope A is a segment and the 
residues at the two endpoints of A arc the classical residues of qy- at zero and 
infinity. We give the general definition in Section 4. 

The coefficients ca are integer numbers which depend only on combinatorics 
of the polytopes Ai, . . . , A„. They can be computed by counting certain complete 
flags of faces of the Minkowski sum A. We discuss this in Section 2. 

A particular case of Theorem 1 produces a formula for the mixed volume of the 
polytopes Ai, . . . , A„ in terms of the vertices of the Minkowski sum A and combi- 
natorial coefficients. Furthermore, Theorem 1 provides with an explicit elimination 
algorithm for systems (1). 

The present paper is expository. Proofs of the results that we mention here are 
contained in [GK02, Kho99, Sop02]. 

2. Combinatorial coefficients 

2.1. Definition. Consider a collection of n convex polytopes Ai, . . . , A„ in R™. 
Let A be their Minkowski sum, A = Ai + • • • + A„. It is not hard to see that each 
face T of A has a unique decomposition into the sum of faces of the Ai: 

F = Fi + ---+r„, a face of A,. 

We will call Tj the i-th summand of T. A face T of A is called locked if at least 
one of its summands has zero dimension, i.e. is a vertex. A vertex A of A is called 
critical if all the proper faces of A that contain A are locked. 

Furthermore, consider a continuous piecewise-linear map x : A — > R" whose 
i-th component X i is non-negative and vanishes precisely on those faces T of A 
whose i-th summand is a vertex. Such map is called characteristic. Note that 
the preimage of the origin under a characteristic map \ consists of the vertices 
of A only. Also in a neighborhood of each critical vertex A the map % sends the 
boundary <9A to the boundary <9R" of the positive octant. 

Definition 1. The combinatorial coefficient ca at a critical vertex A is the 
local degree at A of the restriction of \ to the boundary dA: 

XA :(dA,A)^(dRl,0). 

One can check that it is independent of the choice of a characteristic map \ 
(see [Kho99] for details). 

The sign of the combinatorial coefficient depends on the orientation of R™ and 
the order of the polytopes. Indeed, by definition, the local degree of \A ■ (dA, A) — > 
(<9R™,0) depends on the orientation of the space containing the polytope A and 
on the orientation of the target space. The latter is determined by the order of the 
components Xi of Xi which corresponds to the order of the polytopes. 
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Now we will give the precise definition of generic relative position of n polytopes 
in R n . Let Ai, . . . , A„ be convex polytopes in R™. Every linear functional £ on R™ 
defines a collection of faces Af , . . . , A| of the polytopes such that the restriction of 
£ to Ai attains its minimum precisely at A?. 

Definition 2. We say that n polytopes Ai, . . . , A„ in R™ have generic relative 
position if for any non-zero linear functional £ on K" at least one of the faces 
Af , . . . , A| is a vertex. 

Note that if Ai, . . . , A„ have generic relative position then every proper face 
of the Minkowski sum A is locked. Therefore the combinatorial coefficient ca is 
defined at every vertex A of A. 

2.2. Explicit description. The following theorem gives an explicit descrip- 
tion of the combinatorial coefficient as the number of certain complete flags of faces 
of A, counted with signs. The sign of a complete flag 

r° c r 1 c •• • c c r\ di m r = i 

of faces of a polytope in 1" is the orientation of any frame of vectors (v\, . . . , v n ), 
where Vi starts at T° and points strictly inside P. The sign depends on the orien- 
tation of R™. 

Theorem 2. [Sop02] Let Ai, . . . , A„ be a collection of convex polytopes in R™ 
and A their Minkowski sum. The combinatorial coefficient ca at a critical vertex 
A is equal to the number of all complete flags 

a = r° c r 1 c • • • c i™- 1 c r n = a, dimr = i 

counted with signs, where T l is an i- dimensional face of A whose first i summands 
have positive dimensions and the last n — i summands have zero dimension. 

Notice that this description is highly non-symmetric in the order of the poly- 
topes. In fact, the actual number of complete flags can be different for different 
orderings, whereas the signed number of complete flags is always the same (up to a 
sign), since the combinatorial coefficient is alternating on the order of Ai, . . . , A„ 
as it follows from definition. 

2.3. Example. Two plane polygons have generic relative positions if and only 
if they do not have two parallel sides with the same direction of an inner normal. 
For example, the square and the triangle on Figure 1 have generic relative position. 
Each side of the Minkowski sum comes from either one of the polygons. The 
combinatorial coefficient of a vertex is if both its adjacent sides come from the 
same polygon, and is ±1 otherwise. The sign depends on the orientation defined 
by the two adjacent sides. 

Remark 1. Combinatorial coefficients first appeared in 1996 in [GK96] where 
Theorem 1 was announced. The explicit description (Theorem 2) was found in 
2001 and provided a connection between the Gelfond-Khovanskii residue formula 
and Parshin's reciprocity laws (see [Sop02] ) . It would be interesting to see if this 
description can be useful in mixed volume computation and elimination process 
which we discuss in the next section. 
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Figure 1. The Minkowski sum of two polygons and combinatorial coefficients. 

3. Applications of Theorem 1 

3.1. Mixed Volume. The special case of Theorem 1 when q = 1 gives a 
formula for the number of solutions of the system counting multiplicities. One can 
prove that the residue res^y 3 - A • • • A y*) is independent of the coefficients of the 
polynomials /j and equals det(A l7 . . . , A n ), where Ai is the i-th summand of the 
vertex A (see Section 2.1). On the other hand, according to the BKK theorem the 
number of solutions is nl times the mixed volume of the Newton polytopes. 2 This 
results in a formula for the mixed volume of a collection of n lattice polytopes with 
generic relative position: 

(2) n!V(Ai,...,A„) = (-l) n Y, c A det(A u ...,A n ), 

AGvort(A) 

where the sum is taken over the vertices A of the Minkowski sum A, A, is the i-th 
summand of A, and ca is the combinatorial coefficient at A. 

Let us remark that this formula holds true for non-lattice polytopes as well. 
This was proved by Khovanskii [Kho99] using regular polyhedral subdivisions. 

3.2. Elimination Algorithm. According to classical elimination theory the 
projection of the solution set of the system (1) to the i-th coordinate line can be 
given as the solution set of a single polynomial equation in one variable. In fact, 
Theorem 1 provides us with an explicit algorithm for obtaining such an equation. 

For example, let us obtain a polynomial p{t\) whose roots are the first coor- 
dinates of the solutions of the system. First, compute the number N of solutions 
of the system (counting multiplicities) by taking q = 1. This is the degree of the 
polynomial p. Then taking q = t\, t\, . . . , we get the sums of powers of the 
roots of p. Finally, using the Newton formulas we can express the coefficients of 
the polynomial p via the sums of powers of its roots. To illustrate this procedure 
we give an example in Section 5. 

4. Definition of ves A {q^- A • • • A 

Here we include the definition of the residue res^^y 1 A • 
A of the Minkowski sum of the Newton polytopes of /i , . . . , /, 

2 The Bernstein-Kushnirenko— Khovanskii theorem says that the number of solutions of a 
generic system with fixed Newton polytopes is n! times the mixed volume of the polytopes. 
However, as it was shown by Khovanskii, if the polytopes have generic relative position then 
the number of solutions (counting multiplicities) is the same for all systems with given Newton 
polytopes. 



■Ay*-) at a vertex 
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Before we give the general definition let us look at the one-dimensional case. 
Consider a mcromorphic form q^-, where q, f are Laurent polynomials in one vari- 
able. The possible poles of this form are the non-zero roots of /, zero, and infinity. 
The residue of qj- at a non-zero root a of / is equal to fi(a)q(a), where /i(a) is 
the multiplicity of the root. Since the total sum of residues is zero the sum of the 
values of q over the non-zero roots of / counting multiplicities is the negative sum 
of residues at zero and infinity: 

X! K a )l( a ) = - ( res o <7y + resoo <?y J . 

Recall that reso <7y is the coefficient of 1/t in the Laurent expansion of gy at 

t = and — rcsooq— is the coefficient of 1/t in the Laurent expansion of g— at 

t = oo. We denote them by resA q^f and res_B 3y\ respectively, where [A, B] is the 
Newton polytope (segment) of /. 

Thus we obtain Theorem 1 for n = 1: 



X] M°M°) = - ( res A q y - res s i-j 

Note that the combinatorial coefficients are ca = 1 and cb = — 1- 
The following definition generalizes the notion of the Laurent expansion at zero 
(infinity) to the case of several variables. 

Definition 3. [GK02] Let / be a Laurent polynomial 

/(*) = X \ k t k 

fcez™ 

and A a vertex of the Newton polytope of /. Then Xa ^ and so we can write 
/ = \At A {l — h) for some Laurent polynomial h. Consider a series 

(3) 1 J = J - i L- 1 - r x- A H-* (1 + h + h > + ...,. 

Note that any monomial can appear only in a finite number of terms h l , thus (3) 
is a well-defined power series. Given a Laurent polynomial g, the formal product 
of g and the series (3) is called the Laurent expansion of g/f at the vertex A. 

Now we are ready to define the residue res^gy- A • • • A y 2 -). Note that the 
Newton polytope of the product /i . . . /„ is the Minkowski sum A of the Newton 
polytopes of /i, . . . , f n . Therefore, the Laurent expansion of g / (/i . . . f n ) is defined 
at any vertex A G A for any Laurent polynomial g. 

Definition 4. The number res^qy^- A • • • A y 5 -) is the coefficient of u 1 1 in 
the Laurent expansion of 

<lJf 

fl ■ ■ ■ fn 

at the vertex A. Here J/ is the Jacobian 
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5. Example 

We now consider an example of a system of two polynomial equations in two 
unknowns: 

(4) fi(x,y) = a w x + a 01 y + a 22 x 2 y 2 = 0, f 2 (x, y) = b 00 + b 12 xy 2 + b 21 x 2 y = 0. 

The following are the Newton polygons Ai, A2 of /1, f 2 and their Minkowski 
sum A with the combinatorial coefficients at the vertices. 





By (2) the number of solutions of the given system (4) (which is twice the 
mixed volume of Ai, A2) is equal to 
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6. 



We are going to obtain two degree six polynomials in one variable whose roots 
are the x- and y-coordinates of the solution points of (4), respectively. 

To make our notations shorter we denote the residue res a(q^ A at the 
vertex A = (l,m) by res(; jm ) (q). We compute the residues at the six vertices of A 
for q = x k , 1 < k < 6: 

res ( o,i)(a;' c ) = res (1 , 0) (x fc ) = res (1 , 3) (x fc ) = res( 3;4) (a; fc ) =0, 1 < k < 6 



res (3j i)(x fc ) 



res (4 , 3) (x fe ) 



k = 1,2,4,5 
k = 3 



k 



0, 

36qq (am 621 — ^22 fcoo) 

aio&2i ' 
3b oo ( 2a lQ ^12 (aoi&2l— 2a 2 2 6oo ) + (aoi 621 -ci22boo) 2 ) 

u 10 u 21 

0, 

— 3&i2(aoi&2i— aiofcl2+2a22fcoo) 

022^21 ' 
-3bi 2 ( 2a 22boo(noib2l— 2aiobi2+3a22boo) + (aoifc2l -aiobl2) 2 ) 
n 2 ft 4 ' 

a 22 u 21 



k = 1,2,4,5 
k = 3 

k = 6 



Applying Theorem 1 for q = x k we get the sums Sfc of the fc-th powers of the 
x-coordinates of the solution points of (4): 



si 

S3 



s 2 = s 4 = s 5 = 0, 
resr 



res (3i i)(x 3 ) - res( 4i3) (x 3 ), 
se = rcs (3il) (x 6 ) - res (4i3) (x 6 ). 

We now use the Newton formulas 

s k - sfc^iff! + • • • + (-l)*- 1 *!©*-! + (-l) k kcr k =0, 1 < k < 6 

to obtain the elementary symmetric polynomials a k - 
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0-6 = 



(Ti = (72 = C4 = &5 = 0, 

«3 001^21 (O10&12 + a22&00) - (Ol0&12 - O22&00) 2 

3 aioa 2 202i 
sj - 3s 6 _ a|j6oo6i2 

18 010022^2! ' 

Thus we get a polynomial equation for the x coordinate: 

6 _ aoi&2i(aio&i2 + Q22&00) - (aio&ig - Q22&00) 2 3 Qoifcoo&i2 _ „ 
aioa220 2 i aioa220 2 i 
Since the system (4) is symmetric, an equation for the y coordinate is similar: 

6 _ QlO&12(aQlfr21 + Q22&00) ~ (OQljgl ~ Q^Oo) 2 3 0^00 &21 _ Q 

aoia22&i2 aoia 2 2& 2 2 

Remark 2. The fact that the constant terms of the above polynomials are so 
simple is not a coincidence, but a consequence of a general result due to Khovan- 
skii [Kho99]. Namely, if the Newton polytopes of a system have generic relative 
position then the product of the z-th coordinates of the solution points is always a 
Laurent monomial in the coefficients of the system (see [Kho99] for details). 

Remark 3. Let us remark that the two equations above can be found by 
computing the resultants of fa and fa with respect to x and y, respectively. However 
in higher dimensions computing resultants can be a harder problem than finding 
residues. 
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